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INTRODUCTION 

Certain fiber composites are often modeled as equivalent homogeneous solids [1]. 
Continuum models of this sort are useful in analyzing wave propogation in fiber 
composites when the wavelengths under cons deration are long compared to the mean 
fiber diameter. 

The input-output characterization of a homogeneous transversely isotropic plate is 
investigated by tracing P and SV waves. Following the work in [2], the reflection of a P 
wave at a stress-free plane boundary in a semi-infinite transversely isotropic medium is 
considered first. It is reestablished that an incident P wave reflects a similar P wave and 
an SV wave. It is also reestablished that the angle of reflection is equal to the angle of 
incidence whenever the plane boundary where the reflection occurs is parallel to the 
isotropic plane of the transversely isotropic medium. The angle of reflection of the 
reflected SV wave is found as a function of the angle of incidence of the incident P wave. 

The plane of isotropy of the equivalent transversely isotropic continuum plate lies in 
the midplane of the plate and is parallel to each face of the plate. The P waves experience 
multiple reflections at each face of the plate. At each reflection a P wave and an SV wave 
are reflected back into the medium. The SV waves also experience multiple reflections, 
producing a reflected P wave and a reflected SV wave at each reflection [2]. The 
reflected SV wave is reflected with an angle of reflection equal to the angle of incidence 
of the incident SV wave 12) and the reflected P wave is reflected with an angle of 
reflection equal to the angle of incidence of the incident P wave that produced the 
incident SV wave The amplitude coefficients of the P and SV waves reflected by an 
incident P wave are calculated as functions of the angle of incidence and plotted. Further, 
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the amplitude coefficients of the P and SV waves reflected by an incident SV wave 
produced by mode splitting are calculated as functions of the angle of incidence of the 
incident P wave that produced the incident SV wave and plotted. 

A path notation is then defined to aid in the tracing of series of reflected P and SV 
waves through the medium. Paths are determined by path parameters equal to the number 
of P waves in the path, the number of SV waves in the path and the number of SV wave 
to SV wave reflections in the path. A path amplitude coefficient is defined in terms of the 
path parameters and tabulated. It is found that more than one path may have a given set of 
path parameters and a path multiplicity function is derived to count the number of distinct 
paths with the same path parameters. A net path amplitude coefficient is defined as the 
sum of the amplitudes of all paths with the same total combination of P waves and SV 
waves. The net path amplitude coefficient is tabulated. 

Finally, a theoretical output voltage from the receiving transducer is calculated for a 
tone burst (a periodic input voltage of finite duration) by neglecting the effects of mode 
splitting. 
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REFLECTION OF INCIDENT P WAVE AT STRESS-FREE PLANE 
BOUNDARY IN SEMI-INFINITE TRANSVERSELY ISOTROPIC MEDIUM 


WITH PLANE OF ISOTROPY PARALLEL TO PLANE BOUNDARY 

1. REFLECTED P AND SV WAVES 

A plane progressive stress wave may be represented as 

(u, v, w) = A(P x ,P y , P t ) exp{/co(5 r x + S y y + S t z - r)} ( 1 ) 

for example, see [31, where u , v and w are the displacement components of a point in the 
medium along the x, y and z axes, respectively; A is the amplitude of the particle 
displacement; P x , P, and P, are the components of the unit vector of particle displacement 
along the x, y and z axes, respectively; i « V^T; co denotes radian frequency; S x , S y and S, 
are the components of the slowness vector, which points in the same direction as the 
normal to the wave front and has magnitude equal to the reciprocal of the phase velocity, 
see [4], along the x, y and z axes, respectively; and t denotes time. 

A plane progressive P wave is incident on the plane boundary of a semi-infinte 
linearly elastic transversely isotropic continuum with plane of isotropy parallel to the 
plane boundary. Define a cartesian coordinate system (x, y, z) as follows: the plane 
boundary of the medium contains the x and y axes, and the z axis is the zonal axis of the 
medium (see Fig. 1). The generalized Hooke’s Law, when written relative to the (x, y, z) 
coordinate system, is [4] 
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x„ = C tl ujc + C 12 v ,y + C 13 h\z 
X„ = C n ujc + C u v,y + C, 3 w,z 
X xi — Cj 3 w>x C| 3 v f y *1" (7 33 h*,z 
= C 44 (m,z + w,x) 
x >I = C 44 (v,z+w,>') 

% = C 66 (u,y + vpc) (2) 

where for / = J x tj is a normal stress and for i * j t iy is a shear stress; denotes partial 
differentiation with respect to the succeeding variable; C u , C n , C JV C n and C M are the 
five independent elastic constants for a linearly elastic transversely isotropic medium; 
and is equal to 1/2 (C u - C n ). 

The stresses associated with a plane progressive P wave are evaluated by substituting 
the expression for the displacement components of a point in the medium, given by eqn. 
(1), into eqn. (2) to be 

x„ = /coA (C U SJ> X + C n S y P y + exp{tco(5 jr x + S y y + S,z - r)} 

x„ = UaA {C n Sf x + C u S y P y + C n S,P z ) exp{to)(5^ +S y y + S t z - 1)} 
x„ = i(M (C n S x P x + C n S y P y + C 3J S,P,) txp{i(D(S x x +S y y + S t z - 1)} 
x XI =ioiA (C 44 S,P X + C^SJ 3 ,) exp{/ (DCS,* + S y y + S,z - r)} 

Xy, = iuA(C M S,P y + C u S y P,) exp{i'co(S,x +S y y +S z z -/)} 

x v = icoA (C^ y P x + CJSJ 3 ^ exp{/'co(.S,x + S y y + S,z - r)} (3) 

The stress boundary conditions on a stress free plane boundary require [3] 
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(4) 


t5?+£>=o 


V +<>=(> 

iff + tf-O 


where the Xf> represent the stresses associated with the incident wave and the Xf> 
represent the stresses associated with any reflected waves. The stresses determined by 
eqn. (3) satisfy the boundary conditions, eqn. (4), only when the frequency, 0), of any 
reflected wave is equal to the frequency of the incident wave [3] and 


e(l) = c<*> (5) 

where S (t> and S y !l> ;ire the x and y components of the slowness vector of the incident 
wave, respectively; and S « and S y <*> are the x and y components of the slowness vector of 
any reflected wave, respectively [2], Eqn. (5) establishes that the incident and reflected 
waves lie in the same plane, called the plane of incidence. For computational ease, and 
without loss of generality, assume the plane of incidence coincides with the y-z plane. 
Then the x components of the slowness vectors of the incident and the reflected waves 
vanish, satisfying 

5f = 5f ) = 0 ( 6 ) 

A P wave with slowness surface in a plane containing the zonal axis of a transversely 
isotropic medium possesses a quasi-longitudinal displacement [5]; therefore, a P wave 
traveling in the y-z plane in the (x, y, z) coordinate system has unit vector of particle 

displacement 


(p x ,p,,p,)=(0,pp,p?’ > ) 


(7) 
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where Pp and Pp are, respectively, the directionally dependent y and z components of 
the unit vector of particle displacement of a P wave. 

A wave with P x , P y and P t given by eqn. (7) and S x given by eqn. (6) has stresses 
determined by eqn. (3) satisfying 

v T «;v T «*° (8) 

An SV wave with slowness surface in a plane containing the zonal axis of a 
transversely isotropic medium is quasi-transverse [5]. Therefore, an SV wave traveling in 
the plane x=0 in the (x, y, z) coordinate system has P z , P y and P t according to 

(P x ,P r = (9) 

where P y < sv > and P t (SV) are, respectively, the directionally dependent y and z components of 
the unit vector of particle displacement of an SV wave. 

A wave with unit vector of particle displacment satisfying eqn. (9) and x component 
of slowness vector determined by eqn. (6) has stresses calculated from eqn. (3) satisfying 

^=Ta=0 

v x »;v x »*° (10) 

An SH wave with slowness surface in a plane containing the zonal axis of a 
tranversely isotropic medium possesses a transverse displacement [5]. Thus, for an SH 
wave traveling in the plane x=0 in the (x, y, z) coordinate system the unit vector of 
particle displacement is given by 

(P x ,P y ,P,) = ( 1,0,0) (11) 
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A wave with P t , P } and P, given by eqn. (11 ) and x component of slowness vector 
equal to zero according to eqn. (6) has stresses determined from eqn. (3) satisfying 


x « = x „ = x » = x ,. =0 

x „*0 

T^O < 12 > 

A P wave is incident on the plane boundary of a semi-infinite tranversely isotropic 
medium. The stresses at the point of incidence ihat are associated with the incident P 
wave are constrained by eqn. (8) and must satisfy 

x£> = 0 ( 13 > 

where the superscript ® is used to denote properties associated with the incident wave. 

The stress boundary conditions for any reflected waves are evaluated from eqn. (4) 
and eqn. (13) and can be expressed as 

J» = 0 <14) 

where the superscript <*> is used to denote properties associated with any reflected waves. 
Eqn. (12) indicates that the X a component of the siress tensor associated with an SH wave 
is nonzero. Therefore, no SH wave will be reflected back into the medium by an incident 
P wave, because such a wave would create a nonzero t n stress component and violate the 
stress boundary conditions, eqn. (14). Eqn. (8) indicates that the components of the stress 
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tensor associated with a P wave satisfy the stress boundary conditions and eqn. (10) 
indicates that the components of the stress tensor associated with an SV wave also satisfy 
the stress boundary conditions. Therefore, an incident P wave will result in a reflected P 
wave and a reflected SV wave [2]. 

2. SLOWNESS SURFACE OF P WAVE 

The equations of motion relative to the (x, y, z) coordinate system are [4] 

x a jc + x^y + x„,z = pu,tt 
x xy pc + x yy ,y+x y ,,z=pv,tt 

V+V> + V =pw,« (!5) 

where p is the density of the medium. 

When the components of stress are calculated according to eqn. (3) and the 
displacement components are calculated according to eqn. (1), the equations of motion 
can be written as [2] 

(C U S 2 + CJ>1 + C U S 2 - p )P X + (C u + Cjs x s y p y + (C, 3 + w/, = 0 
(C n + Cjs x s y p x + {CJ>1 + C, x S 2 y + C M S 2 - p )P y + (C, 3 + C M )S y S t P t = 0 
(C„ + C M )S x S t P x + (C t3 + C M )S y S x P y + (C^S 2 + C M S 2 + C n S 2 - p )P t = 0 (16) 

Eqn. (16) can be put into matrix form as 

lB)(P x J> y J>/ = [< 0] (17) 

where [ B ] is a 3x3 matrix with entries 
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( 18 ) 


»ii = (C,A J + C«S, 2 +C« 4 S, 2 - p) 

^22 = "*■ ^1 l^y + ^ 44^1 ~ P) 

6 33 = (C 44 5 i 2 +C 44 5, 2 + C, 3 5 J 2 -p) 

^12 = ^21 = (^12 + CfdSJiy 
b l3 = b il = (C li + C 44 )S x S, 
bn — bj 2 = (C 13 + C^S^S, 

(P J 3 J 3 ) is the unit vector of particle displacement and [0] is the 3x 1 zero matrix. The 
matrix [5], defined by eqn. (18), is a symmetric matrix; therefore, there exist three real 
eigenvalues. One eigenvalue corresponds to an SH wave, another eigenvalue corresponds 
to a P wave, and the third eigenvalue corresponds to an SV wave [2]. 

The plane wave solution is found by setting the determinant of matrix [5], the matrix 
of the coefficients of the unit vector of panicle displacement, equal to zero [6]. 
Expanding the determinant of [fl] and solving for the three roots yields equations for the 
three slowness surfaces. The slowness surface for a P wave is [5] 

i ic, , + c„) (s; + s?) + <c„ + c„)s ; + j [<c, , - c„> (s; + s;> + <c„ - c„>s, 2 f + 

4(S 2 + S 2 )S 2 [(C n - C„) (C„ - CJ - (C„ + C„) ! |}"j = p (19) 

3. ANGLE OF REFLECTION OF REFLECTED P WAVE 

According to eqn. (5) the y components of the slowness vectors of the incident P 
wave, the reflected P wave and the reflected SV wave are equal. Thus eqn. (5) can be 
written as 
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( 20 ) 


b=s^=sp=sf V) 

where b is a constant and S y <'>, Sj F > and S y !SV) are the y components of the slowness vectors 
of the incident P wave, the reflected P wave and the reflected SV wave, respectively. The 
z component of the slowness vector of the reflected P wave, Sj F) , is found from the 
equation for the slowness surface of a P wave, eqn. (19). The values of the x and y 
components of the slowness vector of the reflected P wave, Sj F> and S y < F> , given by eqn. 
(20) and eqn. (5), respectively, are substituted into eqn. (19) and Sj F > is found to satisfy 

5f ) = ±S l (/) (21) 

where Sf> is the z component of the slowness vector of the incident P wave. The slowness 
vector of the incident P wave points out of the medium and the slowness vector of the 
reflected P wave points into the medium (see Fig. 1), therefore the relationship between 
S< F> and S<‘> is 



( 22 ) 


The angle of incidence is defined as the angle between the slowness vector of the 
incident wave and the normal to the boundary at the point of incidence. The angle of 
incidence of the incident P wave, Q ( '>, is then 


0 (/) = arctan 


' S< ;) ' 

-sf 


(23) 


where the point of incidence is taken to be on the plane boundary of the semi-infinite 
body (see Fig. 1). 
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The angle of reflection is defined, in a manner similar to the angle of incidence, as the 
angle between the slowness vector of the reflected wave and the normal to the boundary 
at the point of reflection. The angle of reflection of the reflected P wave, 0"\ is then 


8 (/>) = arc tan 


r*n 



(24) 


Evaluating Q< F > using eqn. (20) and eqn. (22) and comparing the result to eqn. (23) 
establishes that 


0 (/>) = e co (25) 

Thus, the angle of reflection of the reflected P wave is equal to the angle of incidence of 
the incident P wave [2]. 


4. ANGLE OF REFLECTION OF REFLECTED SV WAVE 

The equation for the slowness surface of an SV wave is found from the determinant of 
the matrix [S], defined by eqn. (18), and is [5] 


j(C„ + CJ (S] + S, 2 > + (C„ + C„)S, 2 - f [(C„ - CJ <S 2 + S, 2 ) + (C„ - CJS‘f + 
4(s; + SjtfKC,, - o (C„ - CJ - (C„ + Oll'H = p (26) 


where 5,, S } and 5, are the x, y and z components of the slowness vector of the SV wave, 
respectively; p is the density of the medium and C,., C 1} , C 3J and C u are elastic constants. 
The z component of the slowness vector of the reflected SV wave, S/ iV >, is found from 
eqn. (26) by substituting in the values for S™ and the y and x components of the 



slowness vector of the reflected SV wave, given by eqn. (20) and eqn. (5), respectively. 
The substitution provides two values for Sj sv >, one positive and one negative. The positive 
value corresponds to a wave traveling into the medium and is labeled Sj sv >. The angle of 
reflection of the reflected SV wave, 0 W , is then 


f 


6' 


,(SV). 


arctan 


£(SV)\ 

sT\ 


(27) 


The angle of reflection of the reflected SV wave is, in general, not equal to the angle of 
incidence of the incident P wave [2]. 


A P wave incident, with angle of incidence 0«, in the y-z plane in the (x, y, z) 
coordinate system on the plane boundary of a semi-infinite body has y and z components, 
S y (l > and Sj‘>, of slowness vector satisfying 


cCO 

•(/)_ ? 
tan 0 (/) 


(28) 


Substituting eqn. (5) and eqn. (28) into eqn. (19) allows the equation for the slowness 
surface of a P wave traveling in the plane x=0 to be written as 


A 

2 


(Sf) 2 + 1 (SjY tan" 2 0 (/) + |J[c (Sf) 2 + D (Sf f tan" 2 0 (/) ] 3 


1/2 


4(Sj°) 4 tan -2 0 (/) (CD -E 2 )\ = p 


(29) 


where the conventions: 
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A = C n + C M 
B = C M + C 33 

c = c 1 I -c , 4 

D — C 33 — C M 
E = C, 3 + C 44 

have been used for computational ease. Eqn. (29) can be solved for S y (l> yielding 


( 30 ) 


in 


s?=- 


- + ^ tan- 2 0 (/) + ;{[C+D tan“ 2 0 (, f + 4 tan' 2 0 (;) [CD - E 2 ]} 1 ' 2 


(31) 


An expression for Sj sv > as a function of can be derived from the equation for the 
slowness surface of an SV wave by writing eqn. (26) as 


{ 1 1/2 

Squaring both sides of eqn. (32) and simplifying the resulting expression produces a 
quadratic equation in (S/ SV5 ) 2 . Using the quadratic formula to find the roots of this 
equation produces the expression for (S / SVJ ) 2 



|p S?* + E 2 S? V)1 + |(3pBCD - 2p BE 2 - p AD 2 }S ( y sv9 + 


-3ABCD 


+ ABE 2 + 2C 2 D 2 -3CDE 2 + E" + 


2 . r-4 


b 2 c 2 a 2 d 2 


1 1/2 


\5?* + p 2 D 2 \ 



( 33 ) 
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The z component of the slowness vector of the reflected SV wave is then computed as a 
function of S y <SV) according to 


S™ = jjpfi -^y-Sf v>2 + ^ Sf* + E 2 S^ 9 + |[3p BCD - 2p BE 2 - p AD 2 ]S^ + 

\-3ABi 

L 2 


4 b 2 c 2 a 2 d 2 


+ ABE 2 + 2 C 2 D 2 - 3 CDE* + E* + 


[Sf^ + p 2 d 2 \ 


B^_D^ 
2 2 


, 1/2 


(34) 


Substituting eqn. (20) into eqn. (34) allows Sj sv > to be calculated as a function of the y 
component of the slowness vector of the incident P wave. The value of Sj sv > is then given 
by 


Sf V) = jjpfl -~-sf + ^Y~sf + E 2 Sf v * + j[3pflCD -2pBE 2 -pAD 2 \sf + 


- 3 ABCD 


1/21 


+ ABE 2 + 2 C 2 D 2 - ?>CDE 2 + E* +^-j— sf + p 2 D 2 \ [• + 

4 4 y 


2 2 


, 1/2 


(35) 


Substituting eqn. (20) and eqn. (35) into eqn. (27) allows 0™ the angle of reflection of 
the reflected SV wave to be computed as a function of Sj'> according to 
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6< SV) = arctanj{S,} 




{ [3p BCD - 2pBE 2 - pAD 2 ]Sf ■ 


-1ABCR +AB E 2 + 2C 2 D 2 -3CDE 2 +E* + 

2 4 


B l C l A Z D‘ 
+ — 


\sf + p 2 D 2 


\l/2l 


J J 


B 2 

2 



(36) 


where S « is given as a function of 0 W by eqn. (31 ). 

The values of the material constants of a typical fiberglass epoxy composite are given 
in [4] as C, =10.581x10 N/m 2 , C„=4.679xlO' N/irt*, C n =40.741xl0» N/m*, C 44 =4.422xl0» 
N/m 2 , C m =3.243x 10 9 N/m 2 and p=1850 kg/m 3 . The slowness surface of a P wave traveling 
in the y-z plane in the fiberglass epoxy composite is obtained by substituting the 
numerical values of the constants into eqn. (19) and setting the x component of the 
slowness vector equal to zero. The slowness surface of an SV wave traveling in the y-z 
plane in the fiberglass epoxy composite is found similarly from eqn. (26). The first 
quadrant of the slowness surfaces thus obtained tire shown in Fig. 2. 

The angle of reflection of the reflected SV wave is found as a function of the angle of 
incidence of the incident P wave by substituting the values of the material constants into 
eqn. (36). The angle of reflection of the reflected SV wave is shown in Fig. 3 as a 
function of the angle, of incidence. 
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fiMPiTT.niTTPIJT CHARACTERIZATION OF FIBER COMPOSITE 

1. ANGLES OF INCIDENCE AND REFLECTION OF P AND SV WAVES 

Certain fiber composites may be modeled as a linearly elastic tranversely isotropic 
homogeneous continua [1]. A fiber composite modeled as such a solid in the form of an 
infinite plate where the plane of isotropy lies in the midplane of the plate is to be studied. 
A cartesian coordinate system (x, y, z) is chosen such that the x-y plane coincides with 
the plane of isotropy of the plate and the plate is bounded above by the plane z=h!2 and 
bounded below by the plane z=-hJ2. 

A transmitting transducer and a receiving transducer are assumed to be coupled to the 
top face of the plate in the y-z plane and separated by a distance L. The input electrical 
voltage of the transmitting transducer is a known function of time, V t {t). The output 
electrical voltage of the receiving transducer is an unknown function of time, V (r). The 
transmitting transducer is assumed to be a point transducer that converts an input 
electrical voltage into a stress at a point on the boundary of the plate; the stress then 
travels through the plate as stress waves. The receiving transducer is assumed to be a 
point transducer that converts the stress at a point on the boundary of the plate into an 
output voltage. 

A representative fiberglass epoxy composite plate will be characterised. The compos- 
ite plate has an equivalent continuum model with material properties: 
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p=l 850 kg/m 3 
C, =10.58 lxlO’N/m 2 
C, =4.098x10’ N/m 2 
C, ,=4.679x10* N/m 2 
C 3 3=40.741x10 9 N/m 2 
(^=4.422x1 O’ N/m 2 
C w =3.243xl0 9 N/m 2 

The transducers are assumed to be separated by a distance L=10 cm and the plate is of 
thickness h=5 cm. 

In the following analysis only those stresses at the receiving transducer associated 
with P waves generated by the receiving transducer will be considered. The P waves 
produced by the transmitting transducer experience multiple reflections at each face of 
the plate before reaching the receiving transducer. Since the isotropic plane of the plate 
lies in the midplane of the plate and is parallel to the top and bottom faces of the plate 
each reflection may be treated as the reflection of a plane progressive wave on the plane 
boundary of a semi-infinite linearly elastic transversely isotropic continuum with plane of 
isotropy parallel to the boundary. Therefore, at each reflection the incident P wave 
reflects a P wave and an SV wave; the angle of reflection of the reflected P wave is equal 
to the angle of incidence of the incident P wave [2] and the angle of the reflection of the 
reflected SV wave is given by eqn. (36). 

Assume a P wave is incident, with angle of incidence 0 m , on one of the plane 
boundaries of the plate. An SV wave is reflected back into the medium with angle of 
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reflection 0^. The relationship between the y components, Sf and S y isv >, of the slowness 
vectors of the incident P wave and the reflected SV wave, respectively, is given by eqn. 
( 20 ) as 

5j /) = 5f V) (37) 

The SV wave then travels through the plate and is incident with an angle of incidence 
equal to 0w on the opposite face of the plate (see Fig. 4). An SV wave and a P wave are 
reflected by this new incident SV wave [2]. The reflected SV wave has angle of reflection 
0 (mv> equal to the angle of incidence of the incident SV wave [2], 

q( RSV) = e (W) (38) 

The reflected P wave has angle of reflection 

The relationship between the y components of the slowness vectors of the incident SV 
wave, the reflected SV wave and the reflected P wave is [2] 

S (5l,) = S (RSV) = S (Rf>) (39) 

where Sj sv) , S } < KSV> and Sj* r> are the y components of the slowness vectors of the incident 
SV wave, the reflected SV wave, and the reflected P wave, respectively. Substituting eqn. 
(37) into eqn. (39) demonstrates that 

5«0 _ S (*r) (40) 

The z component, of the slowness vector of the P wave reflected by the incident 
SV wave is found by substituting eqn. (40) into eqn. (1 1) to satisfy 

5 («’ ) = ±5 (/ ) (41) 
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where S » is the z component of the slowness vector of the initial P wave. The P wave 
reflected by the incident SV wave and the initial P wave are both traveling from the upper 
face of the plate towards the lower face of the plate (see Fig. 4), so the relationship 
between S< KF > and S< 1 > is 


S {Rf) = S {I) 

The angle of incidence 0^of the initial P wave is defined by 


(42) 


0 (/> = arctan 


( c<'> A 


-si 


(/) 


(43) 


and the angle of reflection of the P wave reflected by the incident SV wave, is 
defined by 


Q (RP) = arctan 


( \ 

-sr 


(44) 


Substituting eqn. (40) and eqn. (42) into eqn. (44) and comparing the result to eqn. (43) 
establishes that 


0 (^> = 0 (O (45) 

Thus, an incident P wave with angle of incidence 0* ,J will cause a series of SV and P 
waves to be produced by the initial and subsequent reflections. The SV waves will 
always be reflected with angle of reflection 0^, given by equation (36), and the P waves 
will always be reflected with angle of reflection equal to (K 
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2. AMPLITUDE COEFFICIENTS ASSOCIATED WITH INCIDENT P WAVE 


The amplitude coefficients of the P and SV waves reflected by a P wave incident in 
the y-z plane on the plane boundary of semi-infinite transversely isotropic medium are 
given by the pair of simultaneous equations [2] 


and 


sy? + syi n +a (p - p \ spp f + syp) + 

A (P ~ SV) (S {SV) P {SV) + S {SV) P (SV} ) = 0 

+Cnsy* ) +A (p - p xc n syf ) +c^sppy + 
A <r - SV) (c iy sf V) Pf° + C 3 ,Sf>f V) ) = 0 


(46) 


(47) 


where is the amplitude coefficient of the reflected P wave and is the amplitude 
coefficient of the reflected SV wave; S«, Sf> and S™ are the y components of the 
slowness vectors of the incident P wave, reflected P wave and reflected SV wave, 
respectively; S/'\ Sj p > and S< sv > are the z components of the slowness vectors of the 
incident P wave, reflected P wave and reflected SV wave, respectively; PJ'K Pf> and P f <sv > 
are the y components of the unit vectors of particle displacement of the incident P wave, 
reflected P wave and reflected SV wave, respectively; PJ'K Pj F> and P t < sv > are the z 
components of the unit vectors of particle displacement of the incident P wave, reflected 
P wave and reflected SV wave, respectively; and C 13 and C J3 are elastic constants. Eqn. 


(46) can be solved to provide the explicit expression for A <F F> , 




S?p<? + S?pp +A (P ~ SV) {Sf' X 5V) + spm 

SPpP + SPpP 


(48) 


Substituting eqn. (48) into eqn. (47) yields 
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S (l)p{l) + £ (l)p (/) + A (P -SV)^(SV)p(SV) + £(SV)p(SV)^ 

IrFf^fW' x 

[C.tf'P'P+C^f'Pr | ■M'- w tC 1 ^f ) />f' , + CX’O = 0 (49) 

Eqn. (49) is an implicit function for A< psv > dependent only on the material properties and 
the slowness vectors and the particle displacement vectors of the incident and reflected 
waves. This equation can be solved to provide the explicit function for A (PSV > 

A (/> ~ SV) = { - C„syi l) + [(Sf/f + X 

(C^P^+C^SPPf)] + [SPPr + SpP?']} 

{-[(5f>r + (C 13 Sf Vf > + C 33 5f ) Pf>)] + 

[sppp+ S?PP} + C i3 Sf V) Pf V) + C 33 5f (50) 

The expression for A< psv > given by eqn. (50) can be substituted into eqn. (48) to 
produce an expression for A fPP > independent of A< p - sv >. The amplitude coefficient of the 
reflected P wave is then 

A (P ~ P) = -(SpP® + S?p?) (SpPp + spppp - 

<$™P™ + S™PT) (SPpP + Sf >Pf r 1 X 

{-Cvjsfpf - cjtpp? + (sppp + sppP) X 

(Sppp + S^pp + ( C n sppp + c^pp ”)} * 

{-(sr/T + sppn (C l3 sppp+ c 33 spp t (P ) + 

(SPpP + spp t (p ) + C l3 Sf V) Pf V) + C 33 5f >r}"‘ (5 1 ) 

The value of S 3 < r > is found as a function of the angle of incidence 0^ from eqn. (31) to 
be 
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in 


y + f tan" 2 8 (/) +'-{[C+D tan" 2 0 (/) ] 2 + 4 tan" 2 Q l,) [CD - E 2 ] } 1/2 J 

where p is the density of the medium and A, B, C, D and E are functions of the material 
properties defined by eqn. (30). Eqn. (20) establishes the equality of S f (t >, Sj F > and Sj sv > and 
allows S m and S y fsv> also to be found as functions of &'h Substituting eqn. (52) into eqn. 
(28) establishes the value of S<'> as 


S /} = - 


tan 6 (/> 


jj + *tan- 2 e (/) + ;{[C +D tan -2 0 (/> ] 2 + 4 tan " 2 % (, \CD -£ 2 ]} 1/2 J 
Substituting eqn. (52) into eqn. (20) to find S y <sv > and substituting that value into eqn. (32) 
provides the expression for Sj sv> in terms of 0 <l> , as 


tan 0 (Sl ° 


[J+f tan^e^+jdC +£> tan -2 0 (/) ] 2 + 4 tan -2 0 (/) [CD 
where 0w is the angle of reflection of the reflected SV wave and is determined as a 
function of 0 m by eqn. (36). The relationship between S<’> and Sj F > is determined by eqn. 
(22) so the z component of the slowness vector of the reflected P wave can also be found 


as a function of 0^. 
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Thus, the components of the slowness vectors of the incident and reflected waves are 
determined by eqns. (20), (22), (52), (53) and (54) as functions of only the material 
properties and the angle of incidence. 

The y component of the unit vector of particle displacement of a P wave traveling in a 
transversely isotropic continuum is calculated in [5] and can be expressed for a wave 
traveling in the y-z plane as [2] 


P y = 


h p -ds; 


[(//. - DS 2 ) 2 + E 2 S 2 S 2 ] 


(55) 


where the conventions of eqn. (30) have been used. Similarly, the z component of the unit 
vector of particle displacement of a P wave traveling in the y-z plane can be expressed as 


P = 


ESyS, 


[{H p - OS 2 ) 2 + E 2 S 2 S 2 ] ia 


(56) 


The coefficient H p in eqn. (55) and eqn. (56) is defined by [2] 


«,=^cs;+ds, ! +| 


(CS 2 + DS 2 ) 


■4 S 2 S 2 (CD -E 2 ) 


1/2 


(57) 


Substituting eqn. (28) into eqn. (57) allows the coefficient H p to be written as a function 
of only Sf> and the material properties. The H p coefficient associated with the incident P 
wave, Hj'>, is then 
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= tan- 2 0 (/) + 


2 


^C(Sf) 2 + D(Sf) 2 tan" 2 0 (/) j - 


4(5f) 4 tan" 2 e (/) (CD 



Eqn. (58) can be factored as 
where 


( 58 ) 


(59) 


/// = ^| C +D tan~ 2 0 (/) + [(C +£> tan -2 0 (/) ) 2 - 4 tan -2 0 (,) (CD -£ 2 )]' /2 i 


(60) 


Using eqns. (20), (22) and (28), eqn. (57) can be evaluated to provide an expression for 
Hj F >, the H p coefficient associated with the reflected P wave, as 


H p ) = \ j c (Sf ) 2 + D ( 5 f) 2 tan '* 0(0 + (S^ 2 + £> (Sj 0 ) 2 tan -2 0 (/) j 


4(5f > ) 4 tan -2 0 (/) (CD -E 2 ) 


ni/2 


( 61 ) 


Comparing eqn. (61) and eqn. (58) establishes that 

H^ = H ( P (62) 

The y component of the unit vector of particle displacement of the incident P wave is 
computed by substituting eqn. (28) and eqn. (59) into eqn. (55) to be 
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(63) 



(S { ^H p '-D(S^) 2 tan' 2 0 (/) ] + 
[{Sf) 2 H; - D (S { ‘ } ) 2 tan' 2 0 (/> J + 


- 1 - 1/2 


£ 2 (5;Vtan- 2 e (/) 


Canceling a factor of (S/OVCS/9 2 from the right-hand side of eqn. (63) produces an 
expression for P<‘> dependent on 0« and the material properties. The y component of the 
unit vector of panicle displacement of the incident P wave can then be written as 



H/-D tan' 2 0 (/) 

[(/// - D tan' 2 6 (/) ) 2 + E 2 tan' 2 0 (/> ] 1 ' 


(64) 


Similar substitutions allow the z component of the unit vector of particle displacment of 
the incident P wave to be calculated as a function of the angle of incidence and the 
material properties. Eqn. (28) and eqn. (59) are substituted into eqn. (56); a factor of 
(5 m ) 2 /(5/0 2 is canceled from the resulting expression and Pj'> is then determined by 



E ' 
tan 0 (/) - 


_ 1 - 1/2 
(/// - D tan' 2 0 (/) ) + E 2 tan' 2 9 (/) 


(65) 


The components of the unit vector of particle displacement of the reflected P wave are 
found by similar methods. Eqns. (20), (22), (28), (59) and (62) are substituted into eqn. 
(55) to produce an expression for P <*> dependent on only the material properties and 0« 
Comparing the equation for P < F > resulting from the above substitutions to eqn. (65) 
establishes the relationship between P< F > and Pj'> as 

p (P) = pin ( 66 ) 

1 y y 

Substituting eqns. (20), (22), (28), (59) and (62) into eqn. (56) provides an expression 
for P< p > dependent only on Q<’> and the material properties. The value of Pp is given by 
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( 67 ) 


p<n _ . 

*» 


tan 


e (/ >L 


_ " 1-1 (2 

(H/-D tan" 2 0 (/) ) + E 2 tan" 2 0 (/) J 


Comparing eqn. (67) and eqn. (64) establishes 

p (n = _p(n (68) 

The components of the unit vector of particle displacement of an SV wave traveling in 
a transversely isotropic continuum are given in [5] and can be written for the present case 
of a wave traveling in the y-z plane in the (x, y, z) coordinate system as [2] 


Py = 


H^-DS 2 

[(//„ - DS 2 ) 2 + E 2 SyS?] ia 


and 


(69) 


P = 


ES y S, 

[C H„-DS 2 ) 2 +E 2 S 2 y S ?]' 


(70) 


where 


H„ = ^lcS 2 + DS 2 -[(CS 2 ^DS 2 ) 2 -4S 2 S 2 (CD -E 2 )Y 


(71) 


and the conventions defined by eqn. (30) are used. 

Substituting the value of Sj sv> given by eqn. (32) into eqn. (71) allows the H s 
coefficient associated with the reflected SV wave, Hf v \ to be written as 
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1 1 c (>^f V> ) 2 (5f v> ) 2 tan“ 2 0 (SV) - ffc (Sf v> ) 2 -f- Z> (5f V) ) 2 tan -2 0 <SV> Y - 


4(5f v Vtan“ 2 0 (5V) (CD -£ 2 )Jj 


( 72 ) 


where is the angle of reflection of the reflected SV wave and is given as a function of 
the angle of incidence Q"> of the incident P wave by eqn. (36). The right-hand side of eqn. 
(72) can be factored as 


where 



(73) 


H„' = ~ jc + D tan -2 0 <SV) - [(C +D tan 2 0 (W) ) 2 - 4 tan -2 ti (SV) (CD - £ 2 )] 1/2 j (74) 

The y component of the unit vector of particle displacement of the reflected SV wave 
can be calculated by substituting eqn. (32) and eqn. (73) into eqn. (69), then, canceling a 
factor of (S/vy/iS^) 2 from the resulting expression to be 

(SV> //„'-£> tan" 2 0 (5l ° 

P (S } = C75) 

y [(//„' - D tan -2 0 <sl,) ) 2 + £ 2 tan -2 0 (51/) ] 1/2 

Similar substitutions allow the z component of unit vector of particle displacement of 
the reflected SV wave to be written as a function of Q<'> and the material properties. Eqn. 
(32) and eqn. (73) are substituted into eqn. (7()i; a factor of (Sj sv> y/(S y ISV) y is canceled 
from the right-hand side of the ensuing expression and Pj sv > is then given by 
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{HJ - D tan" 2 0 (SV °) 2 + D tan' 2 0 W ' J“ 


(76) 


p(SV) _ . 


tan 0 (SV) 


The amplitude coefficients of the P wave and SV wave reflected by an incident P 
wave can be determined as functions of only the angle of incidence 0 m of the incident P 
wave and the material properties. The value of Sp, the y component of the slowness 
vector of the incident P wave, is found as a function of Q"> from eqn. (52). The y 
components, Sp> and S y <SV) , of the slowness vectors of the reflected P wave and reflected 
SV wave, respectively, are set equal to Sp according to eqn. (20). The value of Sp, the z 
component of the slowness vector of the incident P wave, is determined by eqn. (53) and 
the z component, Sp>, of the slowness vector of the reflected P wave is set equal to the 
negative of Sp according to eqn. (22). The value of S/ SV) , the z component of the slowness 
vector of the reflected SV wave, is found from eqn. (54) where 0 fIV > is the angle of 
reflection of the reflected SV wave and is defined in terms of 0« by eqn. (36). The values 
of Pp and Pp, the y and z components, respectively, of the unit vector of particle 
displacement of the incident P wave, are determined by eqn. (64) and eqn. (65) where the 
coefficient H' is determined by eqn. (60). The y component of the unit vector of particle 
displacement of the reflected P wave, Pp>, is set equal to Pp according to eqn. (66); and 
Pp>, the z component of the unit vector of particle displacement of the reflected P wave, 
is set equal to the negative of Pp according to eqn. (68). The values of P y fSV> and P t <SV) , the 
y and z components of the unit vector of particle displacement of the reflected SV wave, 
are found from eqn. (75) and eqn. (76) where the coefficient is determined by eqn. 
(74). 
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The values of the y and z components of the slowness vectors of the incident P wave 
and the reflected P and SV waves and the y and z components of the unit vectors of 
particle displacement of the incident P wave and the reflected P and SV waves, 
determined as functions of &’>, can then be substituted into eqn. (50) and eqn. (51) to 
obtain an expression for A< p - sv >, the amplitude coefficient of the reflected SV wave, and an 
expression for A (PP >, the amplitude coefficient of the reflected P wave, as functions of Q<'>. 
Thus, the amplitude coefficients of the reflected waves are determined by only the angle 
of incidence of the incident wave and the material properties. 

The amplitude coefficient of the P wave reflected by an incident P wave is calculated 
using the material properties of the representative fiberglass epoxy composite and is 
shown as a function of the angle of incidence of the incident P wave in Fig. 5. The 
amplitude coefficient of the SV wave reflected by an incident P wave is calculated 
similarly and is shown as a function of angle of incidence in Fig. 6. 


3. AMPLITUDE COEFFICIENTS ASSOCIATED WITH INCIDENT SV WAVE 
PRODUCED BY MODE SPLITTING 


The amplitude coefficients and A^v of the SV wave and the P wave, 

respectively, reflected by an incident SV wave are given by the pair of simultaneous 
equations [2] 


£ USV) p ( 1SV ) + £ (ISV)p (ISV) + A ^y-r ytr ) p 


(SV-P) /< ,(RP)r > (RP) n(RP) p (RP) 


+si Kr, pr’)+ 


A (SV - £ (RSV) p (RSV) + £ (RSV) p (RSV}^ _ q 


(77) 


and 
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c 13 sf>r + c n s! m, p, 0S '' > +A ,sv -' > (C„Sl m P;"’ > + CrfTl + 

°+ CjA^r 0 ) = 0 (78) 

where S^' 5 ' 9 and 5/ ,SV) are the y and z components of the slowness vector of the incident 
SV wave, respectively; Sj KSV > and Sj RSV > are the y and z components of the slowness vector 
of the reflected SV wave, respectively; S y <* F > and S,» are the y and z components of the 
slowness vector of the reflected P wave, respectively; P< ISV > and P< rsv > are the y and z 
components of the unit vector of particle displacement of the incident SV wave, 
respectively; P y tRSV) and P, tKSV) are the y and z components of the unit vector of particle 
displacement of the reflected SV wave, respectively; and P y <* F > and P< RPi are the y and z 
components of the unit vector of particle displacement of the reflected P wave, 
respectively. 

Eqn. (77) can be solved for A <SVF> in terms of A< sv -^ as 


A (SV-P) _ VSV)p{ISV)p(ISV) + S (ISV) + A (SV-SV)^SSV)p(RSV) + s (fiSV)p(RSV^ + 


{Sf ) Pf ) + Sf ) P?' ) } 


(79) 


Substituting eqn. (79) into eqn. (78) provides an expression for independent of 

A< sv -r>. Simplifying the resulting equation produces the explicit function for A^^ 

A ( sv - {-C 13 S y SV) P y SV) - C 33 5f>r ) + [(Sf v) pf v) + S ( y SV) Pl ,sv ) x 

(C 13 sf ,) pf ) + CJTPTA + [srrr + > 

{[ -(s^>r + sf ’ /} pd ( c 13 srpr + i *■ 
is;**™ + s^pf^ + [cx sv)p T V) + cjr'pr*) r 1 ( s °) 

Substituting eqn. (80) into eqn. (79) yields an expression for /P svp > independent of 
A< svsv >. The amplitude coefficient of the reflected P wave is 
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Ksrvf 10 + s^r) {c n sf p) pf p) + c^rpr)] [sr p r s r p ^ P) y > x 
{ [-( S, w i», w + 5 f >r v) ) (C^’fVf > + c„sro] + 

[Sf* ’Pf* + S^O + [C I3 Sf V) Pf K) + C 33 sr ) />r > ] }"‘ x 

^(«SK)p(«V) + ^(MV)p(/KV)^ ^(W , ) jD (W , ) + ^(W , )p(R/ , )|- 1 (g^) 

Assume the incident SV wave was produced by the reflection of a P wave incident 
with angle of incidence Q<‘>. Then, the y and z components of the slowness vector satisfy 


S? V) = S? V) (82) 

and 

5 </sv) _ S (SV) ( 83) 

where S y (sv> and S t (SV> are the y and z components, respectively, of the slowness vector of 
the SV wave reflected by the initial P wave. The y component of the slowness vector of 
the SV wave reflected by the initial P wave is determined from eqn. (20) and eqn. (52), 
and Sj svi is determined from eqn. (54). The components of the slowness vector of the 
incident SV wave are thus determined as functions of the angle of incidence of the initial 
P wave, &'>. 

The relationship between the y components of the slowness vectors of the incident SV 
wave and the reflected SV and P waves is given by eqn. (39) as [2] 

S (/5V) = 5 (R£V) = s (RP) (84) 

The z component of the slowness vector of the reflected S V wave is found by substituting 
the expression for S y < KSV > determined by eqn. (84) into the equation for the slowness 
surface of an SV wave, eqn. (26), to satisfy 



(85) 


$(RSV) 


=±s (,SV) 


If the slowness vector of the incident SV wave is pointing in the positive z direction then 
the slowness vector of the reflected SV wave is pointing in the negative z direction and 
vice versa (see Fig. 4); therefore, the relationship between Sj KSV > and Sj lsv > is 


S («SV) = _ S (1SV } ( 86 ) 

where Sj RSV > and Sj ,sv > are the z components of the slowness vectors of the reflected SV 
wave and the incident SV wave, respectively. Substituting eqn. (84) into eqn. (82) and 
then using eqn. (20) and eqn. (52) provides an expression for S y < RSV > dependent on only the 
material properties and 0">. Substituting eqn. (86) into eqn. (83) and substituting the 
resulting expression into eqn. (54) provides an expression for Sj RSV > as a function of the 
material properties and 0w. 

The relationship between S<* r) and the y component, S<'K of the slowness vector of the 
initial incident P wave is given by eqn. (40) as 


Sf ) = Sj /) (87) 

The relationship between Sj* r> and the z component, Sj'>, of the slowness vector of the 
initial incident P wave is given by eqn. (42) as 

( 88 ) 

Substituting eqn. (87) into eqn. (52) produces an expression for S y <* F > dependent only on 
the material properties and &'>, and substituting eqn. (88) into eqn. (52) produces an 
expression for S t mn as a function of 0^. 
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The components of the unit vector of particle displacement of an SV wave traveling in 
a transversely isotropic continuum are given in [5] and are reproduced for the present 
case of a wave traveling in the y-z plane in the (x, y, z) coordinate system in eqn. (69) 
and eqn. (70) following [2]. The incident SV wave is the same wave as the SV wave 
reflected by the initial incident P wave, therefore, the y and z components of the unit 
vector of particle displacement of the incideni SV wave are equal to the y and z 
components, P y ,sv> and Pj sv >, of the SV wave reflected by the incident P wave; 

pilSV) = p(SV) (89) 

and 


p<JSV)_p{SV) (9Q) 

where P y (SV) and Pj sv> are given by eqn. (75) and eqn. (76), respectively. 

The coefficient H n , defined by eqn. (70), can be evaluated for the reflected SV wave 
by substituting eqn. (86) into eqn. (28) to establish 


n(RSV) J y 

1 ~ tan 0 (5l,) 


(91) 


then, substituting eqn. (82) and eqn. (91) into eqn. (71) to evaluate the H n coefficient 
associated with the reflected SV wave, Hj RSV >, as 


H™ = 4 C (Sf v ) 2 + D (Sf vy ) 2 tan- 2 9 W - [jc(Sf V) f + D (Sf V) f tan' 2 0 (S ^ 
4(5f V) ) 4 tan -2 e (SV) (CD -E 2 ^] 


(92) 
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Comparing eqn. (92) and eqn. (74) establishes that 


(93) 

The y component of the unit vector of particle displacement is found by substituting 
eqns. (82), (91) and (93) into eqn. (69) to be 


p(RSV) 


= "(Sf - D tSfV tan" 3 e CJ:VI ] + 

- D (Sf'V tan" 3 0 <SV >J + 4 E\S? V >)‘ tan" 3 ^ 


(94) 


Canceling a factor of (Sj SV) ) 2 /(S y (SV) ) 2 from the right-hand side of eqn. (94) and comparing 
the resulting expression to eqn. (75) establishes that 

p(*sv) _ p(sv ) ( 95 ) 

The z component of the unit vector of particle displacement of the reflected SV wave 
is found by similar substitutions. Eqns. (84), (89) and (93) are substituted into eqn. (70). 
A factor of (S^y/iSj^) 2 is canceled from the right-hand side of the subsequent equation 
to produce an expression for Pj KSV > dependent on only the material properties and 0 WV) , the 
angle of incidence of the incident SV wave. Comparing the expression for P< RSV > thus 
generated to eqn. (76) produces the relationship 

p(RSV) _ _p(SV) ( 96 ) 

The components of the unit vector of particle displacement of a P wave traveling in a 
transversely isotropic contiuum are given in [5] and are reproduced for the present case of 
a wave traveling in the y-z plane in eqn. (55) and eqn. (56) following [2]. The coefficient 
H p , defined by eqn. (57), can be evaluated for the reflected P wave by substituting eqn. 
(88) into eqn. (28) producing the expression 
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(97) 


_c<0 

c(W) 

' tan0 (/) 


then, substituting eqn. (87) and eqn. (97) into eqn. (57) to evaluate the H p coefficient 
associated with the reflected P wave, H< RF >, as 


= ^jcCS^ 2 + D (Sf) 2 tan' 2 0 (/) + ^CtfJ 0 ) 2 + D tan" 2 0 

-II/2] 

4(5; /) ) 4 tan' 2 0 </) (Cf> -£ 2 )J 


(98) 


Comparing eqn. (88) and eqn. (60) establishes that 

(99) 

The y component of the unit vector of particle displacement of the reflected P wave is 
found by substituting eqns. (88), (97) and (99) into eqn. (55) to be 


p(RP) _ 


(S^H/ - D (Sf) 2 tan -2 0 (/) ] + 

((S«)V - D (S { ;) 2 tan“ 2 0 (/) J + E 7 (S^) 2 tan" 2 0 (,) 


(100) 


Canceling a factor of (S/ SV} )V(S/ V 0 2 from the right-hand side of eqn. (100) and comparing 
the resulting expression to eqn. (64) produces the relationship 

p(RP) _ p (i) (101) 

y y 
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The z component of the unit vector of particle displacement of the reflected P wave is 
found by similar substitutions. Eqns. (88), (97) and (99) are substituted into eqn. (56). A 
factor of {Sf v >) 2 /(S y < sv >) 2 is canceled from the right-hand side of the subsequent expression 
for Pj F > and the resulting expression is compared to eqn. (65) to establish that 

piftr) _ p(D (102) 

The amplitude coefficients of the SV wave and the P wave reflected by the incident 
SV wave produced by mode splitting can then be determined as functions of only the 
angle of incidence of the incident P wave, that produced the incident SV wave, and 
the material properties. The value of S y <lsv >, the y component of the slowness vector of the 
incident SV wave, is found as a function of 0w by substituting eqn. (82) into eqn. (52). 
The y components, Sj RSV) and Sj* F> , of the slowness vectors of the reflected SV wave and 
reflected P wave, respectively, are set equal to Sj RSV> according to eqn. (84). The value of 
Sj lsv >, the z component of the slowness vector of the incident SV wave is found by 
substituting eqn. (83) into eqn. (54), where 0^ is defined as a function of B (,) by eqn. 
(36). The z component, Sf™, of the slowness vector of the reflected SV wave is set equal 
to the negative of SJ ISV> according to eqn. (86). The value of Sj* F> , the z component of the 
slowness vector of the reflected P wave, is found by substituting eqn. (88) into eqn. (53). 
The values of Pj ISV) and P< tsv >, the y and z components of the unit vector of particle 
displacement of the incident SV wave, are found by substituting eqn. (89) into eqn. (75) 
and substituting eqn. (90) into eqn. (76), respectively. The coefficient Hj is determined 
as a function of &'> by eqn. (74). The y component, P } < KSV >, of the unit vector of particle 
displacement of the reflected SV wave is found by substituting eqn. (95) into eqn. (75), 
and P< KSV >, the z component of the unit vector of particle displacement of the reflected SV 
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wave, is found by substituting eqn. (95) into eqn. (76). The y and z components, Pw and 
P<*n, of the unit vector of particle displacement of the reflected P wave are found by 
substituting eqn. (101) into eqn. (64) and substituting eqn. (102) into eqn. (65), 
respectively. 

The values of the y and z components of the slowness vectors of the incident SV wave 
and the reflected SV and P waves and the y and z components of the unit vectors of 
particle displacement of the incident SV wave and the reflected SV and P waves, 
determined as functions of Q ( '>, can then be substituted into eqn. (80) and eqn. (81) to 
determine A WVSV) , the amplitude coefficient of the reflected SV wave, and A< SVF >, the 
amplitude coefficient of the reflected P wave, as functions of (R Thus, the amplitude 
coefficients of the waves reflected by an SV wave produced by mode splitting are 
determined in terms of only the angle of incidence of the P wave that produced the 
incident SV wave and the material properties. 

The amplitude coefficient of the P wave reflected by an incident SV wave produced 
by mode splitting is computed using the material properties of the representative 
fiberglass epoxy composite. The amplitude coefficient of the reflected P wave is shown 
as a function of the angle of incidence of the initial P wave in Fig. 7. The amplitude 
coefficient of the SV wave reflected by an incident SV wave produced by mode splitting 
is calculated similarly. The amplitude coefficient of the reflected SV wave is shown as a 
function of the angle of incidence of the initial P wave in Fig. 8. 
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4. WAVE PATHS 


It is assumed that the receiving transducer does not produce an output voltage unless 
there exists a series of reflected waves beginning (in this case) with a P wave produced 
by the transmitting transducer, and ending with a P wave incident at the receiving 
transducer. Let such a series of reflected waves be called a path from the transmitting 
transducer to the receiving transducer. Both P waves and SV waves experience mode 
splitting when incident on the stress-free plane boundary of a transversely isotropic 
continuum [2]; therefore, a path may include both P and SV waves. 

A simple path from the transmitting transducer to the receiving transducer may be 
represented as 

M = P P SV P (103) 

Eqn. (103) should be read as follows: M is the path from the transmitting transducer to 
the receiving transducer that begins with a P wave produced by the transmitting 
transducer that is incident on the lower face of the plate, producing a second P wave that 
is, in turn, incident on the upper face of the plate, producing an SV wave that then is 
incident on the lower face of the plate, producing a final P wave that is incident at the 
receiving transducer. The complexity of this description demonstrates the usefulness of 
the path notation. 

A more complicated path of demonstrative value is 

N = P P SV P SV SV P P (104) 
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Path N is interesting because it includes all four possible types of reflection: P wave 
producing P wave, P wave producing SV wave. SV wave producing P wave, and SV 
wave producing SV wave. These four types of reflection have amplitude coefficients 
A (F P> , A (FSV> , A ,SVF) and A (SV ~ SV) . 

5. NUMBER OF P AND SV WAVES IN PATH: WAVE INDICES 

When the initial P wave of a path from the transmitting tranducer to the receiving 
transducer is incident with a known angle of incidence, 0">, on the lower face of the plate 
the angles of incidence and reflection for all subsequent P waves and SV waves in the 
path are also known. The angle of reflection and the angle of incidence of any P wave in 
the path are both equal to 0«, and the angle of reflection and the angle of incidence of any 
SV wave in the path are both equal to 0™, which is defined as a function of 0» by eqn. 

(36). 

A path from the transmitting transducer to the receiving transducer is therefore 
characterized by the angle with which the initial P wave of the path is incident on the 
lower face of the plate. Because the plate is of thickness h, the total distance r p traveled 
between reflections by any P wave in the path is (see Fig. 9) 

_ h (105) 

r '~cose<'> 

The distance b p traveled in the y direction between reflections by the P wave is 

b p = r p sin 0 (/t < 106) 

Using eqn. (105) to evaluate r f in eqn. (106), the distance b p can be written as 
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(107) 


b p —h tan 0 (/) 

The distance r n traveled between reflections by any SV wave in the path is (see Fig. 9) 


h 

~ cos 0 (SV) 


(108) 


where 0w is defined as a function of 0» by eqn. (36). The distance b n traveled in the y 
direction between reflections by the SV wave is 

b„ = r„ sin 0 (JV) (109) 

Evaluating r n in eqn. (109) from eqn. (108) allows the distance b n to be written as 

b JV = h tan 0 (SV) (HO) 

The transmitting transducer and the receiving transducer are assumed to lie in the 
plane x=0 and be separated by a distance L in the y direction; therefore, the total distance 
traveled in the path in the y direction must be L. Thus, if there are r P waves in the path 
and s SV waves in the path, then the wave indices r and s must satisfy the relation 

rbp + sb^^L (1H) 

Substituting eqn. (107) and eqn. (110) into eqn. (Ill) and dividing both sides of the 
resulting expression by h allows the relationship between r and s to be written as 

r tan 0 (/) + s tan S* 510 = ^ 012) 

h 


Since the first wave in the path is produced by the transmitting transducer and every 
other wave in the path is produced by a reflection in the path the total number of 
reflections in the path must be one less than the total number of waves in the path, r + s. 
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Furthermore, the transmitting transducer and the receiving transducer are assumed to be 
coupled to the same face of the plate; therefore, the total number of reflections in any 
path from the transmitting transducer to the receiving transducer must be an odd number. 
Thus, the total number of waves in the path must be an even number, mathematically 
expressed as 

{1,2,3,...} (113) 

where e denotes "is a member of. 

In order for a path to be of interest, the characteristic angle of the path, 0®, the number 
of P waves in the path, r, and the number of SV waves in the path, s, must satisfy eqn. 
(112) and eqn. (113). 

6. PATH AMPLITUDE COEFFICIENT AND REFLECTION INDEX 

When the characteristic angle of a path, 0"', is known, the four amplitude coefficients 
associated with the path, A lFP >, A' sv ^ and A^*\ can be calculated from eqns. (50), 
(51), (80) and (81), respectively. A path amplitude coefficient. A, can then be defined as 
the ratio between the amplitudes of the final P wave and the initial P wave of the path, for 
a nondispersive and nonattenuating medium. This path amplitude coefficient can be 
calculated as a function of 0® according to 

a = {a^~ sv y (A *- f y t a {sv - sv) t (A (sv y o 14) 

where / is the number of P wave to SV wave reflections in the path, m is the number of P 
wave to P wave reflections in the path, n is the number of SV wave to SV wave 
reflections in the path, and p is the number of SV wave to P wave reflections in the path. 
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When the path is written as a list of P’s and SV’s following the notation of eqn. (103), / is 
the number of times the character string P SV occurs in the list, m is the number of times 
the character string P P occurs in the list, n is the number of times the character string 
SV SV occurs in the list, and p is the number of times the character string SV P occurs in 
the list. The path amplitude coefficient for the path defined as N, A(N), is then found 
from eqn. (104) and eqn. (114) to be 

A(N) = (A (p - SV} f(A (p - p} ) 2 (A (SV - s y (A {SV - p ) 2 (1 15) 

The characteristic angle, 0®, together with the wave indices r and s do not uniquely 
determine the path amplitude coefficient. More than one path may have the same 0«, r 
and s; furthermore, paths with the same 0®, r and s may or may not have the same 
amplitude coefficients. The three paths P p P 2 and P 3 defined by 

P 3 =P SV P SV P P 
P 2 = P SV SV P P P 

P, = PPP SVSVP (116) 

have the same wave indices, r=4 and 5=2, and therefore, the same characteristic angle, 
calculated from eqn. (112), but are distinct paths. The path amplitude coefficients, A(P 2 ) 
and A(Pj), of paths P 2 and P 3 , respectively, are equal and computed as 

a(p,) = a(p 3 ) = (A ( '- sv y y F -y (A < sv - s y ( A (SV -y a n) 

but path P, has a different path amplitude coefficient, A(P,), computed as 

a(p 1 )=(A« , - s yy p -y(A (sv - n ) 2 (ns) 

Therefore, another index, a reflection index n, must be found to uniquely determine the 
path amplitude coefficient. 
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The path amplitude coefficient of any path is determined by the characteristic angle of 
the path and the number of each type of reflection in the path. The last wave in any path 
of interest is a P wave; therefore, every SV wave in the path must reflect in the path. The 
number of SV wave to SV wave reflections in the path summed with the number of SV 
wave to P wave reflections in the path must therefore be equal to the total number of SV 
waves in the path, 5 . If there are n SV wave to S V wave reflections in the path, then there 
must be s-n SV wave to P wave reflections in the path. Since the first wave in the path is 
a P wave and the last wave in the path is a P wave, the number of P wave to SV wave 
reflections in the path must be equal to the number of SV wave to P wave reflections in 
the path, s-n. The total number of reflections in the path is equal to r+s-1, one less than 
the total number of waves in the path. Since the number of SV wave reflections in the 
path is equal to the number of SV waves in the path, s, the number of P wave reflections 
in the path must be equal to r-1, one less than the number of P waves in the path. 
Therefore, the number of P wave to P wave reflections in the path must be equal to 
(r-l)-(s-n). The path amplitude coefficient A w can then be calculated as a function of the 
characteristic angle 0 W , the wave indices r and s, equal to the number of P and SV waves 
in the path, respectively, and the reflection index n , equal to the number of SV wave to 
SV wave reflections in the path, according to 

7. MULTIPLICITY FUNCTION 

The number of distinct paths from the transmitting transducer to the receiving 
transducer that have characteristic angle 0 m , wave indices r and s, and reflection index n, 
is found by applying elementary counting theory to the lists of the characters P and SV 
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that represent the paths. The lists are partitioned into character strings, where a character 
string is an ordered subset of a list, of the form P SV...SV and character strings of the 
single character P. The character strings of the form P SV...SV are chosen such that the 
final SV in the character string is followed by a P in the list. The single character 
character strings are chosen such that the P is followed by a second P in the list. Thus, 
each character string in a partition begins with a P and ends with the character preceding 
the next P in the list. Therefore, if a P is followed by a P, in the list, the first P becomes a 
single character character string in the partition; otherwise, a character string consisting 
of the first P followed by as many SV’s as occur in the list, until the next P is 
encountered, is included in the partition (see Table 3 for examples). 

Let the number of character strings, in a list, of the form P SV... SV, containing exactly 
<7+1 SV’s, be represented by k f . Thus, is the number of character strings in the path 
representing one P wave to SV wave reflection and q SV wave to SV wave reflections. 
The number of character strings in a list of the form P SV is then given by k 0 , since in this 
case <7=0. Let the number of character strings of the single character P be denoted by k p . 
The number of SV wave to SV wave reflections is exactly n, so k f must be zero for all q 
greater than n. Furthermore, the total number of SV wave to SV wave reflections in the 
path is equal to 

IV? 

<7 * 1 

so the multiplicity coefficients k t must satisfy the relation 

X k q = n (120) 

, = i 
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The number of SV waves in the path is equal to 


+ 1 ) 

< 7=0 

and must be exactly s. Therefore, the k q s must also satisfy 

i* ? (<7 + l) = 5 (121) 

< 7*0 

In order for the receiving transducer to be excited, the last wave in the path must be a 
P wave. Therefore, k p must be at least one; mathematically, this can be represented by 

k p Z 1 (122) 

The number of P waves in the path is equal to 

k„+ I k 

< 7=0 

and is exactly r. Therefore, the multiplicity coefficients, k q and k p , satisfy 

k p +ik=r (123) 

<?=o 

The number of distinct paths from the transmitting transducer to the receiving 
transducer ending in a P wave, that have r P waves, s SV waves, n SV wave to SV wave 
reflections, and multiplicity coefficients k p , k^ v k 0 and k p , satisfying eqns. (120), (121), 
(122) and (123), is equal the number of distinct lists with k m objects of type 1, k^ x objects 
of type 2, ..., k 0 objects of type n+1 and ^-1 objects of type n+2. The value of k p is 
reduced by one to remove one object of type n+2 from random placement in the list in 
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order to place it at the end of the list; guaranteeing the last wave in the path is a P wave. 
The number v of the distinct lists with k, objects of type i for each of n+2 types of objects 
is given by [7] 


v = 


fi + 2 


X kj 


i * 1 


k x , &2, * • • y k n + j , k n +2 


where 


/ 


m 


i * » 1 % a - 1 * 

is the multinomial coefficient defined by 


(124) 


m 


‘ • •» ^a- 1> 


0 if n { < 0 for any i 

if n, > 0 for all i 


ml 


(125) 


[n 1 !n 2 !.../i a _ 1 !/i a ! 

where x! represents the factorial expansion of x, (x)(x-l)(x-2)...(2)(l). 
The number N of distinct paths is then 


N = 


( " N 

(*,-l) + Ik, 

< 7=0 


(126) 


V 


1 > ^ 0 * • * • ’ - 1 > ^» 

Substituting eqn. (120) and eqn. (123) into eqn. (126) allows the expression for N to be 
simplified to 


N = 


r- 1 

k p ~ ko, • . M k n _i, k t 


(127) 
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The total number of paths, with r P waves, s SV waves and n SV wave to SV wave 
reflections, is found by summing eqn. (127) over all possible values of the multiplicity 
coefficients. Using the constraints of eqns. (120), (121), (122) and (123) allows the total 
number of distinct paths N tMt to be calculated as 


N. 


TJA 


1 C 

: 1 
k m = 0 k a 


«2 

I 

*2 = 0 ! 


f 


r- 1 


k p -hk ( 


0’ 


Ic Ic 


(128) 


where 


a ,= 


n - X k,i 

i = q + 1 


(129) 


where LcJ denotes the greatest integer less than x. 



The characteristic angle, Q<'>, the number of P waves, r, the number of SV waves, s , 
and the number of SV wave to SV wave reflections, n, have been held constant 
throughout the counting process; therefore, all paths counted by eqn. (128) have the same 
path amplitude coefficient. 

The total number of distinct paths with wave indices r and s and reflection index n, 
N and the associated path amplitude coefficient, A w , are tabulated in Table 1 for all 
paths with ten or fewer transits of the plate for the fiberglass epoxy composite. The 
characteristic angle is found from eqn. (112) by fixing r and s and substituting in the 
values of L and h, L= 10 cm and h=5 cm. The values of A w are found from eqn. (1 19) and 
are observed from Table 1 to be either greater than or less than zero for different values 
of r, s and n. This is due to the fact that A***-**, the ratio of the amplitude of particle 
displacement of the SV wave reflected by an incident SV wave to the amplitude of 
particle displacement of the incident SV wave, is greater than zero and A<*-*> and 
A< SVI ‘> are all less than zero (see Figs. 5,6,7 and 8). A negative amplitude coefficient 
indicates that the direction of the particle displacement due to the reflected wave, at the 
point of reflection, is in the opposite sense of the direction of the particle displacement 
due to the incident wave. Whereas, a positive amplitude coefficient indicates that the 
direction of the particle displacement due to the reflected wave, at the point of reflection, 
is in the same sense as the direction of the particle displacement due to the incident wave. 

8. PHASE VELOCITY AND TIME DELAY 

The distance r traveled between reflections by any P wave in a path with 
characteristic angle 0« is given eqn. in (105). If there are exactly r P waves in the path 
then the total distance R p traveled by P waves in the path is 
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(133) 


R ' 


= r 


h 

cos 0 (,) 


where h is the thickness of the plate. 

The distance r n traveled between reflections by any SV wave in a path with 
characteristic angle Q<” is given by eqn. (108). If there are exactly s SV waves in the path 
then the total distance R„ traveled by SV waves in the path is 


R " 5 cos6 (SV) 


(134) 


where & SV) is defined as a function of by eqn. (36). 

The directionally dependent phase velocity of a P wave traveling in a transversely 
isotropic continuum, C p (0), is, for example, see [4] 


f 


C,(0) = 


C44 + C u sin 2 0 + C 33 cos 2 0 + VcT) 

2p 


1/2 


where 


(135) 


G = [(C, , - CJ sin 2 0 + (C M - C 33 ) cos 2 0] 2 + 4(C 13 + C„f sin 2 0 cos 2 0 (136) 

The directionally dependent phase velocity of an SV wave traveling in a transversely 
isotropic contiuum, C^(0), is, for example, see [4] 


Cc V (0) = 


f C44 + Cj , sin 2 0 + C33 cos 2 0 - VGY" 


V 


2p 


where G is defined by eqn. (136). 


(137) 
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The time delay t rJ is defined as the time taken for the final P wave of a path with r P 
waves and s SV waves to reach the receiving transducer after the path is initiated by the 
transmitting transducer. The time delay is then 

'- = c^ + cj^ (138) 

Thus, the time delay is a function of only the geometry of the transducer arrangement, the 
characteristic angle, 0^, the wave indices, r and s, and the material properties of the plate. 

9. NET PATH AMPLITUDE COEFFICIENT 

Eqn. (138) establishes that all paths from the transmitting transducer to the receiving 
transducer containing r P waves and s SV waves experience the same time delay in 
reaching the receiving transducer. Furthermore, all paths with r P waves and s SV waves 
will be in phase at the receiving transducer. Therfore, the amplitudes of all paths with 
wave indices r and s may be added algebraically at the receiving transducer. A net path 
amplitude coefficient A rj may be defined as the ratio between sum of the amplitudes of 
the final P waves of all paths with r P waves and s SV waves and the amplitude of the 
initial P wave of any path with wave indices r and s, when the medium through which the 
waves are traveling is nondispersive and nonattenuating. 

The net path amplitude coefficient is found by summing the product N rjA A rjjt over all 
possible values of the reflection index n, where N rjji is the total number of paths with r P 
waves, s SV waves and n SV wave to SV wave reflections and is given by eqn. (128); 
and A w is the path amplitude coefficient of any path with wave indices r and s and 
reflection index n and is given by eqn. (119). The bounds on n are found by considering 
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the exponents of the amplitude coefficients in eqn. (119). Clearly, if s is equal to zero 
there can be no SV wave to SV wave reflections, and for all values of s greater than zero 
there must be at least one P wave to SV wave reflection. Therefore, n must satisfy the 
relationship 


n € < 


{0} if 5=0 


(139) 


[{0,1,. ..,5 - 1} if s > 0 

where e denotes "is a member of. The lower bound on n is computed by noting that the 
exponents of A <FF > and A ISV - SV > must be nonnegative for all values of n. The lower bound on 
n is then 


n > max (5 + 1 -r,0) 


(140) 


The net path amplitude coefficient is then 


(Ar p r x ^ 5=0 




. £ a h„Af*» if * >0 

=* nwx(j + 1 — rj 0) 


(141) 


The net path amplitude coefficient A r0 = {A < FF) ) rl corresponds to the path P P P...P 
containing (r-1) P wave to P wave reflections and no SV waves. 


The net path amplitude coefficient A rj is shown in Table 2 for all paths with ten or 
fewer transits of the plate for the fiberglass epoxy composite. The value of A rj is found by 
summing the product A rj J^! tjji over all values of n. The maximum and minimum values of 
n are found from eqn (139) and eqn. (140) and the values of A and N are taken from 
Table 1. 
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10. ASSUMPTIONS ON TRANSDUCERS 


The transmitting transducer and the receiving transducer are assumed to be longitudi- 
nal transducers that transform an electrical voltage into a uniform longitudinal stress or a 
uniform longitudinal stress into an electrical voltage. The following approach parallels 
that of [8], 

If an input voltage V i of amplitude V and frequency co is applied according to 

V i = V exp{-i'cor} (142) 

where and t denotes time, the stress t that is introduced into the medium at the 

transducer-medium interface by the transmitting transducer is 

Ty(r) = F x {(a)V exp{-t (cot + <)>,)} (143) 

where F,( to) is the frequency dependent transduction ratio for the transmitting transducer 
in transforming a voltage into a stress and <}>, is a phase angle. In eqn. (142) and eqn. 
(143) the harmonic character of the signals is expressed in complex notation, but only the 
real parts of these and subsequent equations should be considered. The amplitude T of the 
applied stress is then 

T=F l (to)V (144) 

Similarly, if a stress wave producing a stress component xj of amplitude T and 
frequency co that impinges on the receiving transducer is defined as 

x t '{t) = rexp{-tCDf } (145) 

then the output voltage V from the receiving transducer is 

V a (t) = / r 2 (a>)7 , 'exp{-/(G)f + <J» 2 )} (146) 
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where F 2 (w) is the frequency dependent transduction ratio for the receiving transducer in 
transforming a shear stress to a voltage, and <j> 2 is a phase angle. Thus, the amplitude V' of 
the output electrical voltage is 

V' = F 2 ((0)T' (147) 

The characteristics of F t ( co) and F 2 (to) are unknown except that the product /\(a>)F 2 (o>) is 
dimensionless. 

11. DIRECTIVITY FUNCTIONS 

The directivity functions of the stresses associated with P and SV waves traveling in 
semi-infinite transversely isotropic continua are evaluated from the far-field asymptotic 
solutions of the displacement components produced by a harmonic point load imbedded 
in an infinite body in [2]. The directivity functions D*f> of the T. stress associated with P 
waves whose slowness vectors are confined to the y-z plane and are produced by a point 
load located on the plane boundary of a semi-infinite body and acting in the K direction 
in the (x, y, z) coordinate system are reproduced. 

Define a coefficient X R according to 

_ //,S 2 + //,S y 2 + //,S 2 (148) 

I K J 

where I I denotes the "magnitude of’ 

K h = ■ HjS y S y - H,SX) + 

(2 HJS, • HJ y (HJ x S, ■ HJS y S,-H,S x S y ■ H,S,S t )] (149) 

and 
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•■f 


^ S 2 + ^. (S 2 + S 2 ) _ l 


y(5 z 2 + 5;) + y5 f 2 


C44 + C 


13 


\ 2 


Sf(Sx+Sy) 


(150) 


Z denotes the sum with respect to cyclic permutation of S x , S y and S t ; denotes 
differentiation with respect to the succeeding variable carried out assuming S x , S y and 5, 
are independent variables; C u , C 13 , C 33 and are elastic constants of the material; and p 
is the density of the material. 

The directivity function D r J F > of the T„ shear stress associated with a P wave whose 
slowness vector is confined to the plane x=0 and that is produced by a harmonic point 
load located on the plane boundary of an infinite half-space and acting in the y direction 
in the (x, y, z) coordinate system is 



| c 3 3 5 > 2 ^ 

2n Ip * P 


(151) 


where (0,S\S/) is the point on the slowness surface of a P wave where the normal to the 
slowness surface is parallel to the line connecting the point load to the point at which the 
directivity function is to be evaluated, (0 denotes the radian frequency of the harmonic 
point load and X M is evaluated at (OJ/f;). The directivity function D r f> of the x a normal 
stress is 


„<,) WXXfCnC^, C33C* 4 r . 2 „ ' 
=_ ~2k p J " P ' 13 


(152) 
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The directivity function D z y j p > of the x n shear stress associated with a P wave produced 
by a harmonic point load acting in the z direction in the (x, y, z) coordinate system is 


z y ) _ ^(0^;A>C 4 ^ y Y C l , o | C33-C 44 -C 1 3 5 > 2_^ 

Dyt 2n [p y p ’ , 


(153) 


Similarly, the directivity function D z a < F) of the t m normal stress is 


of"- 


X n (0,5 y yS, )tiiS, (C n C i3 Cii ^33„*2 


2k 


Cl 


\ 


+yS, -C 33 


(154) 


J 


The directivity functions D K f> of the stresses associated with a P wave whose 
slowness surface is confined to the plane x=0 and that is produced by a harmonic point 
load can be evaluated at a point M in a semi-infinite body by finding the point (0,5/A/) 
on the slowness surface of a P wave where the normal to the slowness surface is parallel 
to the line TJM connecting the point load to the point M (see Fig. 10). The slope m of the 
line TJM is found from Fig. 10 to be 


Ay 1 

m = — = -77, 

Az tan 0 (/) 


(155) 


The normal to the slowness surface of a P wave is found at any point in the y-z plane 
by setting the x component of the slowness vec tor, S x , equal to zero and differentiating 
the equation for the slowness surface, eqn. (19), with respect to 5 f , the z component of the 
slowness vector. Solving the subsequent equation for the slope, -dSJdS t , of the normal 
yields 
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jr r 

-T - = ) BS, + [QCD - 2 E 2 )S 2 S, + D 2 S, 3 ] x 

do a 

[C 2 sj + (6 CD - 4E 2 )S 2 S 2 + D 2 S*] V1 • + 

^AS y + [C 2 S] + (3 CD - 2 E 2 )S y S 2 ] x 

[C 2 S 4 y + (6CD - 4E 2 )S 2 S 2 + D 2 S*]~ m i (156) 

where the conventions of eqn. (30) have been used. Setting eqn. (155) equal to eqn. (156) 
establishes the constraint on 5/ and 5/, the point on the slowness surface where the 
normal is parallel to the line connecting the transmitting transducer to the point M, 

-^=\bs; + [.icd-2e')s-X + d's-:]* 

[cV + (6 CD - 4E ! )S; 2 S; 2 + D’S," J j * 

Ls; +[c 2 s; ! + <3 cp - 2 £ ! >s;s; 2 ] x 

[c’s; + (6 CD - 4E‘)SX + D X* }'] ( 1 57) 

The point (0 JS y ‘£‘) is then found by setting S x equal to zero in eqn. (19) and solving eqn. 
(19) and eqn. (157) simultaneously. 
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12. STRESS FIELD RADIATED BY TRANSMITTING TRANSDUCER 

If there were no bottom boundary to the plate and the stress waves were propogating 
in an infinite half-space, the path from the transmitting transducer to the receiving 
transducer with characteristic angle 0^/' consisting of 2m P waves and no SV waves 
would travel to a point M^ 0 , in a time r^ 0 , traveling a distance R p2m (see Fig. 11). The 
amplitude of the hypothetical stress x. at the point M^ 0 is and is defined as [2] 


T m 

1 2m ,0 


D 

- t ~r 


k(p) r 
1>2m exd-o/? 


P 2m 


l 


P7m\ 


(158) 


where T is the magnitude of the longitudinal stress generated by the transmitting 
transducer in the z direction, is the value of the appropriate directivity function 

evaluated for the point M^ 0 and a is the P wave attenuation constant of the medium. The 
P wave path, however, does not propagate in an infinite half-space, but instead 
experiences 2m- 1 reflections. The amplitude T„ im 0 of the stress at the receiving transducer 
is thus obtained by modifying eqn. (158) as 


TD*z&rr-' 


°2m,0 


R 


P2m 


expj-aR 


P 7m 


(159) 


where A^ p - f) is the P wave to P wave reflection coefficient, given by eqn. (51), evaluated 
for a characteristic angle of 0*, o m . 

The amplitude t' O2m0 of the output voltage from the receiving transducer is obtained 
from eqn. (147) and eqn. (159) and is 




°2m,0 


R 


P2m 



(160) 
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Substituting eqn. (144) into eqn. (160) yields 


’ °2m.O — R 

^Plm 

where V is the amplitude of the input voltage. Introducing a possible electrical signal 
amplification factor, K, eqn. (161) can be written as 


exp-j-atf 


P2m 


(161) 


V = 

°2m,0 


KF l (<D)F 2 (w)VD*2( A w ) ) 2m 


expS 


-a R 


P2m I 


x P2m 


( 162 ) 


The amplitude T 0tji of the stress at the receiving transducer associated with a path with 
characteristic angle Q r J !> consisting of r P waves and s SV waves can be understood 
qualitatively by visualizing the waves as traveling in a hypothetical multi-layered 
half-space [8] (see Fig. 12). It is assumed that each layer is identical to the original plate 
and that the layers are bonded together in a manner such that an incident wave in one 
layer produces no reflected waves, but produces transmitted P and SV waves in the next 
layer and that the transmission coefficients are the same as the reflection coefficients in 
the original plate. The amplitude of the stress is then set equal to Tj*> the value of the 
stress amplitude at {joint M rj in the hypothetical half-space. The value of T r J u> is then 


T m _ 

1 Tfi 


rP>,A, 

Rp r + R„, 



(163) 


where T is the magnitude of the stress generated by the transmitting transducer in the z 
direction, 15* rj is an equivalent directivity function for a path with wave indices r and s, 
A rj is the path amplitude coefficient, given by eqn. (119), R Pr is the total distance traveled 
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by P waves in the path, R„ t is the total distance traveled by SV waves in the path, a is the 
P wave attenuation constant of the medium and p is the SV wave attenuation constant of 
the medium. 

If is much greater than A rj for all s * 0 and the assumptions 

D k ~D k(P) (164) 

>J V 

and 

a~P (165) 

are then made, it is apparent that the stresses at the receiving transducer associated with 
paths containing P and SV waves are much smaller than the stresses at the receiving 
transducer associated with paths containing only P waves. Therefore, in the input-output 
characterization of fiber composites which may be modeled as transversely isotropic 
contimuum plates, the effects of mode splitting may be neglected. Only those paths from 
the transmitting transducer to the receiving transducer consisting of 2m P waves and no 
SV waves will produce significant stresses at the receiving transducer when eqn. (164) 
and eqn. (165) are satisfied. 

13. OUTPUT VOLTAGE DUE TO TONE BURST 

Assume the input voltage V. to the transmitting transducer is a periodic function of 
time, of center frequency to and duration t r Mathematically the input voltage can be 
expressed as the sum of two periodic functions of frequency to that are 180° out of phase, 
one beginning at time r=0 and one beginning at time t=t r The input voltage is then given 
by 
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V, = V'exp{/G)r}U(f)-l / exp{/G)} U(f - 0 


(166) 


where V is the amplitude of the input voltage, i-V^T , and U(jc) is the unit step function 
defined as 

0 if x < 0 

U(jc) = | (167) 

1 if x > 0 

The periodic nature of the signal in eqn. (166) and subsequent equations is expressed in 
complex notation, but only the real part of the signal is to be considered. 

The voltage output V a at the receiving transducer is found by considering the two 
terms on the right-hand side of eqn. (166) independently and superposing the voltage 
output associated with each input signal. Consider the input voltage V' given by the first 
term on the right-hand side of eqn. (166) 

V' = V exp{/(or} U(r) (168) 

The stress waves produced by a transmitting transducer with input voltage V' will have a 
periodic character with respect to time. Therefore, the amplitude of the output voltage 
associated with each wave path will also be periodic in time and the total output voltage 
V' must be found by the superposition of the contributions of each wave path. 

Only those paths with 2m P waves and no SV waves will be considered. Eqn. (112) 
can then be rewritten as 

2/ntane£ = | (169) 
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where 0^® is the characteristic angle of the paih, L is the separation of the transducers 
and h is the thickness of the plate. Eqn. (169) indicates that the characteristic angle of the 
plate must satisfy 


f 

0^ = arctan 


L 

2mh 


\ 

J 


(170) 


All the trignometric functions of 0^® can be evaluated from eqn. (170) (see Fig. 13). 

A P wave path with 2m P waves will be in phase at the receiving transducer with a P 
wave traveling in a semi-infinite half-space to a point M m (see Fig. 14). The T stress 
associated with a plane progressive P wave is given by eqn. (3) and can be calculated at 
point M m as 




M, 


exr 


CO 


S y L - S. 2mh - 1 

J m 


(171) 


where T Mm is the amplitude of the stress and S ym and S, m are the y and z components of the 
slowness vector of the P wave, respectively. The stress x tjo at the receiving transducer is 
then 



S y L-S, 2 mh - 1 

y*i l m 


(172) 


where T„ m is the amplitude of the stress and is given by eqn. (159). The output voltage V 0m ' 
at the receiving transducer is found by substituting eqn. (172) into eqn. (146) to be 



where <}> is a phase angle. The output voltage is then calculated using eqn. (144) and eqn. 
(159) and including a possible electrical signal amplification factor K as 


V ' = 

°m 


KF l (u)FmVD*2(A£<~ P> ) 2m " ' 




P2m 


expj-aR 


P2m 


X 



(174) 


where V is the amplitude of the input voltage, D*,/’*. is the appropriate directivity 
function, is the P wave to P wave reflection coefficient and R p2m is the total distance 
traveled by the wave path. The total output voltage V' is found by summing eqn. (174) 
over all m and retarding each contribution by the delay time to be 


V/= I 

m = 1 


- RF l (o))F 2 ((o)V D^(A^ n ~ P) ) 2m 


R 


ex P)-° t/? P2*r x 


P 2m 


l 




u (r-U 


(175) 


where t ^ is given by eqn. (138) with s= 0. 

The output voltage V” associated with the input voltage characterized by the second 
term on the right-hand side of eqn. (166) is found similarly to be 


„ - -KF^)F 2 «o) VD^ iA^r-' f t 

v"= X * exp -ctfjx 

m = l K p2m [ 


exp|i|^a)^5 >)(i L - S lm 2mh - 1 j + <}>jj U(r - (/, + r^)) 


(176) 


where t i is the duration of the input signal. 
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The output voltage V B associated with the input voltage V. determined by eqn. (166) is 
found by adding eqn. (175) and eqn. (176) to be 

- DS2(ASr n f m -' f „ 1 

v, = KF,(oi)F 2 «o)VZ i — cxp>, -aR p J/y 


exps i 


(0 




2 mh -rj + <|>]j {U(r - Lj - U(r - (t, + hJ) } (177) 


where D K f>^ is given by eqn. (151), (152), (153) or (154), AJ** is given by eqn. (51), 
R p7m is given by eqn. (133), S ym is given by eqn. (31), S, m is given by eqn. (28), t ^ is given 
by eqn. (138), t i is the duration of the input signal and the characteristic angle Qj'> is 
given by eqn. (170). The sum in eqn. (177) does not include the effects of mode splitting, 
therefore the predicted output voltage V will not exactly equal the true output voltage. 
However, the effects of mode splitting have been shown to be negligible for the initial 
wave paths. Therefore, the initial predicted output voltage should correspond closely to 
experimental data. 


14. SUMMARY OF IMPORTANT RESULTS 

The angle of reflection of an SV wave reflected by an incident P wave is determined 
as a function of the angle of incidence of the incident P wave by eqn. (36). 

The angle of reflection of a P wave reflectttd by an SV wave produced by mode 
splitting is found to be equal to the angle of incidence of the incident P wave that 
produced the incident SV wave, eqn. (45). 

The amplitude coefficients of the P and SV waves reflected by an incident P wave, 
Aff-f) and A< F - SV >, are found as functions of the angle of incidence of the incident P wave by 
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eqn. (50) and eqn. (51), respectively. 

The amplitude coefficients of the P and SV waves reflected by an incident SV wave 
produced by mode splitting, A< SVF > and A< sv ^ v >, are found as functions of the angle of 
incidence of the incident P wave that produced the incident SV wave by eqn. (81) and 
eqn. (80), respectively. 

A path amplitude coefficient, A fJA , is defined for any path from the transmitting 
transducer to the receiving transducer in terms of wave indices r and s, defined as the 
number of P waves in the path and the number of SV waves in the path, respectively, and 
a reflection index n, defined as the number of SV wave to SV wave reflections in the 
path, eqn. (114). 

The number of paths having the same wave and reflection indices, N , are counted in 
accordance with eqn. (127). 

A net path amplitude coefficient, A rj , equal to the product A rj N rj * summed over all 
values of the reflection index for a given pair of wave indices is computed in eqn. (141). 
The product A rjj N rjji can be physically interpreted as the sum of the amplitudes of all 
paths having wave indices r and s and reflection index n. 

The magnitude of the longitudinal stress at a point in a hypothetical multilayered 
half-space due to all paths having the same wave indices is determined by eqn. (163) for 
all time after the arrival of the first path having those wave indices. 

An output voltage from the receiving transducer is determined as a function of time by 
eqn. (177). 
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CONCLUSION AND DISCUSSION 

The input-output characterization of certain fiber composite plates can be studied 

using a transversely isotropic continuum model when the wavelengths under considera- 

♦ 

tion are long compared to the mean fiber diameter. For wavelengths close to, or less than, 
the mean fiber diameter, a continuum plate model is no longer appropriate and the 
inhomogenaities of the composite must be considered. 

In the analysis of a transversely isotropic continuum plate an incident P wave is found 
to reflect a P wave and an SV wave during each reflection at the top or bottom face of the 
plate. The angle of reflection of the reflected P wave is equal to the angle of incidence of 
the incident P wave and the angle of reflection of the reflected SV wave is uniquely 
determined by the angle of incidence of the incident P wave. When the reflected SV wave 
is in turn incident on the opposite face of the plate, a P wave and an SV wave are again 
reflected. This second reflected SV wave is reflected with an angle of reflection equal to 
the angle of incidence of the incident SV wave and this second reflected P wave is 
reflected with an angle of reflection equal to the angle of incidence of the original 
incident P wave. Thus, tracing a series of reflects P and SV waves through the medium 
is simplified because all the reflected P waves art; reflected with angle of reflection equal 
to the angle of incidence of the initial P wave and all reflected SV waves are reflected 
with angle of reflection equal to the angle of reflection of the initial reflected SV wave. 
This simple relationship exists because the isotropic plane of the plate is parallel to the 
top and bottom faces of the plate. If this special geometry is not present a more 
complicated analysis considering each reflection individually is needed to trace a series 
of reflected waves through the plate. 
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A path from the transmitting transducer to the receiving transducer may contain both P 
and SV waves (no SH waves are produced during the reflections). The path amplitude 
coefficient is found from the product of the reflection coefficients encountered at each 
reflection. Because more than one path may have the same number of P and SV waves, a 
net path amplitude coefficient is found by summing the path amplitude coefficients of all 
paths having the same total combination of P waves and SV waves. Table 2 shows the 
magnitudes of the net path amplitude coefficients associated with paths containing no SV 
waves to be an order of magnitude or more larger than the magnitudes of the net path 
amplitude coefficients associated with paths containing SV waves. 

Examining Figs. 5, 6, 7 and 8 reveals that the magnitudes of the reflection coefficients 
for P wave to P wave and SV wave to SV wave reflections are generally close to one, 
especially for the small angles of incidence common for paths from the transmitting 
transducer to the receiving transducer. On the other hand, the magnitudes of the reflection 
coefficients for P wave to SV wave and SV wave to P wave reflections are much smaller 
than one for small angles of incidence. Therefore, paths containing both P waves and SV 
waves, which thus must contain at least one P wave to SV wave reflection and one SV 
wave to P wave reflection in order to start and end with P waves, should be expected to 
have path amplitude coefficients with magnitudes much less than the magnitudes of the 
path amplitude coefficients associated with paths containing no SV waves. Table 1 
reveals that for those values of the wave indices and reflection index for which the 
mutiplicity function is large the magnitude of the path amplitude coefficients tend to be 
much much less than one. Therefore, the increased multiplicity of paths containing both P 
waves and SV waves does not cause a significant increase in the magnitude of the net 
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path amplitude coefficient. 

The magnitude of the net path amplitude coefficients associated with paths containing 
P and SV waves can be expected to be much smaller than the magnitude of the path 
amplitude coefficients associated with paths containing only P waves. Therefore, the 
input-output characterization of a transversely isotropic continuum plate by P waves can 
be carried out neglecting the effects of mode splitting and tracing only those paths 
containing no SV waves. 

This study enhances the theoretical understanding of the nondestructive evaluation 
(NDE) of transversely isotropic media such as certain fiber composites. It also provides 
impetus for further study on the input-output characterization of transversely isotropic 
media by P waves by simplifying the analysis. 
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Table 2 Net path amplitude coefficient A r g as function 
number of P waves r and number of SV waves s for 






Table 3 Partitions, wave indices, reflection index and 
multiplicity coefficients for sample wave paths 


Wave Path 

Partition 

r s n kp k Q ki k 2 | 

P P SV SV P SV P P 
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P SV SV 
P SV 
P 
P 

mm 

P SV SV SV P P P P 

P SV SV SV 
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p 

p 
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P P SV P SV P SV P 

p 

P SV 
P SV 
P SV 
p 

5 3 0 2 3 0 0 

P SV P P SV SV P P 

P SV 
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P SV SV 
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5 3 13 110 
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Fig. 1 Angle of incidence of incident P wave 9^ \ 
angle of reflection of reflected P wave 
9* P) and angle of reflection of reflected 
SV wave 6 shown in y-z plane of (x,y/Z) 
coordinate system. 
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Angle of Reflection of SV wave (degrees) 




Angle of Incidence of P wave (degrees) 


Fig. 3 Angle of reflection of reflected SV wave as f 

angle of incidence of incident P wave for P w 
on stress-free plane boundary of continuum pi 





Fig. 4 Incident P wave produces SV wave which 
experiences inode splitting at opposite 
face of plate. 
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Angle of Incidence of P wave (degrees) 

Amplitude coefficient of SV wave reflected by incident P 
wave as function of angle of incidence of incident P wave 
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Angle of Incidence of P wave (degrees) 

Amplitude coefficient of SV wave reflected by incident SV 
wave produced during mode splitting as function of angle 
of incidence of initial P wave. 



Fig . 7 



Angle of Incidence of P wave (degrees) 

Fig. 8 Amplitude coefficient of P wave reflected by incident SV 

wave produced during mode splitting as function of angle 
of incidence of initial P wave. 
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S V . 

Transmitting , , . S , T > Receiving 

Transducer H *1 0 Transducer 



h A 
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Fig. 9 Wave path p. SV- SV. . . P* SV« P from transmitting 
transducer to receiving transducer containing 
both P and SV waves. 
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Point Load 



P wave traveling to point M in semi-infinite 
body . 






Fig. 12 Schematics of (a) wave path P • SV- SV- P • SV- P 

in continuum plate and (b) wave path 
p. SV* SV-P • SV- P visualized in hypothetical 
multilayered half-space. 
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